flD-fll65  916 
UNCLASSIFIED 


CALIBRATED  L-BAND  TERRAIN  HEASUREHENTS  AND  ANALYSIS  1/2 

PROGRAH  (U>  ENVIRONHENTAL  RESEARCH  INST  OF  HICHIGAN 
ANN  ARBOR  R  U  LARSON  ET  AL  FEB  86  ERIH-1742ee-5-P 
RA0C-TR-8S-266  F19628-84-C-0881  F/G  8/6  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAI  BUREAU  OF  STANDARDS  1963  A 


DIB  fiLc  CQBl  AD-A165  9 


RADC-TR.85-266 
Interim  Report 
F»brucNy  1986 


CAUmm  L-BAKO  TERRAIK 

mimiiEims  m  amalysis  msm 


Environmental  Research  Institute  of  Michigan 


R.  W.  Larson 
E.  S.  Kaslschko 
A.  Moffott 


DTIC 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTIOR  URLIMITED 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griff iss  Air  Force  Base,  NY  13441-5700 


^4 


I 

4 

i 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
I  it  will  be  releasable  to  the  general  public.  Including  foreign  nations. 

I  RADC-TR-85-266  has  been  reviewed  and  is  approved  for  publication. 

I 


APPROVED:  |(  V>/ 

K.  V.  N.  RAO 
Project  Engineer 


APPROVED: 


y 

ALLAN  C.  SCHELL 

Chief,  Electromagnetic  Sciences  Division 


FOR  THE  COMMANDER: 


JOHN  A.  RITZ 

Plans  &  Programs  Division 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  ple?»se 
notify  RADC  (EECT)  Hanscom  AFB  MA  01731-5000.  This  will  assist  us  in  main¬ 
taining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


UNCLASSIFIED 

SeCURITV  CLASSIflCATIONOFT^PAGE 


AMUsn 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 
UNCLASSIFIED 

2a  SECURI’Y  classification  AUTHORITY 

N/A _ 

2b  declassification /DOWNGRADING  SCHEDULE 

n/a 

4  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

174200-5-P 


6a  NAME  OF  PERFORMING  ORGANIZATION 
Environmental  Research 

Institute  of  Michigan 

6c  ADDRESS  (City,  State,  and  ZIP  Code) 

P.O.  Box  8618 

Ann  Arbor  MI  48107 


6b  OFFICE  SYMBOL 
(If  applicable) 


1b  RESTRICTIVE  MARKINGS 
N/A 

3  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

RADC-TR--85-266 

7a  NAME  OF  MONITORING  ORGANIZATION 

Rome  Air  Development  Center  (EECT) 

7b  ADDRESS  (Ofy,  Stafe,  and  ZIP  Code) 

Hanscom  AFB  MA  01731-5700 


8a  NAME  OF  FUNDING /SPONSORING 
1  ORGANIZATION 

Rome  Air  Development  Center 
8c  ADDRESS  forv.  State,  and  ZIP  Code) 

Hanscom  AFB  MA  01731-5000 


8b  OFFICE  SYMBOL 
(If  applicable) 

EECT 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F19628-84-C-0081 


10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  PROJECT 

ELEMENT  NO  NO 

62702F  4600 


PROJECT 

TASK 

NO 

NO 

4600 

15 

11  TITLE  (Include  Security  Clatsification) 

CALIBRATED  L-BAND  TERRAIN  MEASUREMENTS  AND  ANALYSIS  PROGRAM 


WORK  UNIT 
ACCESSION  NO 

73 


12  PERSONAL  AUTHOR(S) 

R.  W.  Larson,  E.  S.  Kasischke,  A.  Maffett 


lia  TYPE  OF  report 
Interim 

16  supplementary  notation 
N/A 


13b  TIME  COVERED 
FROM  Apr  84  TO  As 


14  DATE  OF  REPORT  (Year,  Month.  Day) 
_ February  1986 


PAGE  COUNT 


1 7 _ _ COSATI  CODES _  18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

field _ GROUP _ SUBGROUP  Electromagnetic  Scattering^  ^Ground  Truth;  / 

17 _ 02 _ Terrain  Cross  Section  Measurement^  J  Soil  Moisture  <  •k^L.  - 

20  14  Dielectric  Constant  ^  -  '  (See  Reverse')  / 

19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  j 

This  report  presents  the  initial  results  of  a  three  year  program  for  the  collection  and 
analysis  of  L-band  calibrated  backscatter  data  from  various  terrain  types.  During  the  first 
year  advantage  was  taken  of  two  SAR  data  collection  programs  to  obtain  ground  measurements 
simultaneously  with  the  SAR  flights.  These  data  collection  flights  were  conducted  in  North 
Carolina  during  April  1984  sponsored  by  USGS  and  during  August  -  September  1984  in  Long 
Island,  sponsored  by  ONR.  Ground  data  were  collected  at  specific  sites  during  both  SAR 
flights  and  Included  measurement  of  dielectric  constant,  surface  roughness  and  moisture. 

The  reduction  of  ground  data  and  the  calibration  of  the  SAR  are  described.  Calibration  of 
the  SAR  has  been  verified  and  an  error  of  less  than  3  dB  in  absolute  calibration  of  the 
L-band  channel  demonstrated. 


20  QISTRIBUT'ON/ AVAILABILITY  OF  ABSTRACT 

□  'jnclassified/unlimited  si  same  as  RPT 

22*  NAME  OF  RESPONSIBLE  INDIVIDUAL 
K.V.N.  Rao 


21  ABSTRACT  SECURITY  CLASSIFICATION 
□  otic  USERS  UNCLASSIFIED 

22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 
(617)  861-3106  RADC  (EECT) 

SECURITY  CLASSIFICATION  OF  i-HIS  PAGE 

UNCLASSIFIED 


DO  FORM  1473, 84  MAR 


B3  APR  edition  may  be  uted  until  exhausted 
All  other  editions  are  obsolete 


UNCLASSIFIED 


17.  COSATI  CODES  (Continued) 
Field  Group 

09  03 


18.  SUBJECT  TERMS  (Continued) 

L-band 
SAR  Data 


UNCLASSIFIED 


PREFACE 


i 


The  work  described  in  this  report  was  conducted  by  personnel  in 
the  Radar  Division  of  the  Environmental  Research  Institute  of 
Michigan.  The  work  reported  covers  the  period  April  1984  through 
April  1985  and  discusses  the  results  of  the  first  year  of  a  three- 
year  program.  This  work  is  supported  by  the  Electronic  Systems 
Division,  Air  Force  Systems  Command,  USAF,  Hanscom  AFB,  MA  under 
Contract  No.  F 19628-84-0081 .  The  Air  Force  technical  monitor  is  Dr. 
K.V.N.  Rao,  Hanscom  AFB.  The  SAR  data  utilized  in  the  research  dis¬ 
cussed  herein  were  collected  for  programs  sponsored  by  the  U.S. 
Geological  Survey  (USGS)  under  Contract  No.  14-08-0001-21748  and  the 
Office  of  Naval  Research  (ONR)  under  Contract  No.  N00014-81-C-0692. 

The  principal  investigator  for  this  work  is  Richard  W.  Larson. 
The  work  described  in  this  report  was  performed  by  R.W.  Larson,  A. 
Maffett  and  E.S.  Kasischke. 
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CALIBRATED  L-BAND  TERRAIN  MEASUREMENTS  AND  ANALYSIS  PROGRAM 

1 

INTRODUCTION 


This  report  presents  the  initial  results  of  a  three-year  program 
being  conducted  by  the  Environmental  Research  Institute  of  Michigan 
(ERIM)  for  the  Electronics  Systems  Division  of  the  U.S.  Air  Force 
Systems  Command.  The  goal  of  this  program  is  to  collect,  process 
and  statistically  analyze  calibrated  L-band  clutter  data. 

This  report  gives  the  results  from  the  first  year  of  this  three- 
year  program.  It  also  outlines  the  procedures  which  will  be  used  to 
further  reduce  and  analyze  the  calibrated  L-band  terrain  data.  This 
program  was  designed  to  use  existing  calibrated  synthetic  aperture 
radar  (SAR)  data,  as  well  as  SAR  data  collected  during  two  flight 
programs  in  1984.  During  these  programs,  ground  measurements  of 
surface  roughness  and  dielectric  characteristics  were  obtained  simul¬ 
taneously  with  the  SAR  overflights.  SAR  data  from  one  of  these  data 
collection  programs  has  been  digitally  processed  into  imagery,  and 
using  an  existing  calibration  algorithm,  L-band  terrain  clutter  sta¬ 
tistics  generated.  Ground  measurements  made  in  test  areas  during 
the  data  gathering  period  have  been  reduced  to  obtain  surface  rough¬ 
ness  statistics  and  dielectric  constant  values. 


This  report  is  divided  into  six  chapters,  including  this  intro¬ 
duction.  Chapter  2  presents  an  executive  summary,  which  includes  a 
discussion  of  the  overall  approach,  the  specific  tasks  of  the  three- 
year  program,  and  the  significant  achievements  of  the  past  year. 
Chapter  3  discusses  the  sources  of  data  used  for  this  program. 
Chapter  4  presents  the  data  analysis  and  reduction  techniques  used 
to  produce  calibrated  L-band  terrain  measurements  from  L-band  SAR 
imagery  and  outlines  the  techniques  used  to  reduce  the  ground  mea¬ 
surements.  Chapter  5  presents  the  data  analysis  and  reduction  tech¬ 
niques  to  be  used  during  the  remainder  of  this  program.  Finally, 


Chapter  6  discusses  recommended  research  areas  for  analysis  of  L-band 
terrain  measurements  beyond  the  end  of  this  program.  An  appendix  to 
this  report  presents  the  results  from  the  ground  measurements,  in¬ 
cluding  data  products  from  reduction  of  the  surface  roughness  mea¬ 
surements  and  photographs  of  test  sites. 


Approach  to  the  Calibrated  L-Band  Terrain  Measurements 
and  Analysis  Program 


2 

EXECUTIVE  SUMMARY 


This  report  discusses  the  results  of  the  first  of  a  three-year 
program  to  collect,  process  and  analyze  L-band  clutter  data.  The 
goals  of  this  program  are  to: 

-  Determine  the  statistical  distribution  of  L-band  clutter  data 
as  a  function  of  terrain  types, 

-  Evaluate  existing  L-band  clutter  or  scattering  models,  and 

-  Correlate  L-band  clutter  measurements  with  terrain  scattering 
descriptors. 

The  program  was  designed  to  address  the  above  objectives  using 
existing  or  recently  collected  calibrated  SAR  imagery.  Figure  1 
outlines  the  overall  approach  of  the  program. 

The  program  contains  six  distinct  tasks: 

1.  SAR  data  collection  and  ground  measurements, 

2.  SAR  data  processing, 

3.  Ground  data  reduction, 

4.  Determination  of  the  statistical  distribution  of  L-band 
clutter  data, 

5.  Evaluation  of  L-band  clutter/scattering  models,  and 

6.  Correlation  of  L-band  clutter  data  to  ground  measurements. 

Table  1  summarizes  the  elements  of  each  of  these  tasks  and  pre¬ 
sents  a  schedule  for  their  completion. 

During  the  past  year,  advantage  was  taken  of  two  SAR  data  collec¬ 
tion  missions  which  were  conducted  using  the  CV-580  SAR  system. 
These  missions  were  carried  out  near  Durham,  North  Carolina  during 
April  1984  and  at  and  around  the  Peconic  River  Airport  on  Long 


Island,  NY  during  late  August  and  early  September  1984.  The  North 
Carolina  mission  was  sponsored  by  the  U.S.  Geological  Survey  (USGS), 
and  the  Long  Island  mission  was  sponsored  by  the  Office  of  Naval 
Research  (ONR).  Extensive  SAR  calibration  and  ground  measurements 
were  obtained  coincident  with  these  data  collection  missions.  These 
latter  efforts  were  supported  by  this  program  and  are  described  in 
this  report. 

Descriptions  of  the  SAR  calibration  reference  reflector  arrays 
are  given  along  with  examples  of  calibration  relationships  for  the 
SAR  L-band  channel.  As  many  as  21  reference  reflectors  were  used  to 
calculate  each  SAR  calibration  curve  for  L-band.  Results  from  the 
initial  calibration  verification  analysis  indicate  that  a  calibration 
accuracy  of  less  than  3  dB  was  achieved.  Initial  results  obtained 
from  applying  calibration  data  to  determine  scattering  coefficients 
are  given.  Average  values  of  calculated  (HH  polarization,  63* 
incidence  angle)  included  0.005  for  grass  fields,  0.015  for  forested 
sites  and  0.020  for  corn  fields. 

Ground  measurements  obtained  included  surface  height  profiles, 
dielectric  constant  values,  moisture  content  and  plant  density. 
Surface  height  profile  data  have  been  used  to  obtain  height  spectra, 
slope,  radii  of  curvature  statistics  and  correlation  length.  Exam¬ 
ples  of  these  data  are  presented. 

An  example  of  the  statistical  analysis  procedure  is  given  along 
with  preliminary  results.  An  inverse  Gaussian  (IG)  distribution  was 
selected  to  model  the  histogram  of  the  scattering  coefficients.  A 
three  parameter  description  has  been  derived  for  an  IG  distribution 
using  emperical  data. 

Future  plans  include  continuation  of  the  analysis  using  the  I.G. 
distribution  to  determine  distinct  signatures  that  can  be  related  to 
the  surface  measurements.  Also,  surface  profile  and  electrical  prop¬ 
erties  of  the  surface  will  be  used  to  evaluate  models  that  describe 
the  scattering  characteristics  observed  from  SAR  data. 


Several  surface  scattering  models  have  been  selected  from  the 
literature  as  candidates  for  initial  comparison  to  measured  scatter¬ 
ing  data.  Each  model  utilizes  a  particular  description  of  the  scat¬ 
tering  surface  that  can  be  derived  from  available  measured  ground 
data.  Future  analysis  to  verify  or  to  justify  modifications  to 
existing  models  will  be  conducted.  SAR  data  from  selected  areas  will 
be  used  to  extend  the  range  of  incident  angles  over  which  the  scat¬ 
tering  response  is  available. 


The  approach  taken  to  achieve  the  objectives  of  this  program  was 
to  obtain  L-band  radar  cross  section  statistics  using  data  collected 
by  the  ERIM/Canadian  Centre  for  Remote  Sensing  (CCRS)  CV-580  SAR 
System.  During  the  past  year,  terrain  roughness  characteristics  were 
obtained  simultaneously  with  two  SAR  data  collection  missions  being 
performed  for  the  U.S.  Geological  Survey  (USGS)  and  the  Office  of 
Naval  Research  (ONR).  In  addition,  archived,  calibrated  L-band  SAR 
data  are  also  available  for  determining  clutter  statistics  of  terrain 
types  not  present  in  the  above  data  sets. 

In  this  chapter,  we  will  first  discuss  the  characteristics  of 
the  CV-580  SAR  System  and  the  steps  taken  to  achieve  calibration  for 
this  system.  This  will  be  followed  by  a  description  of  the  two  test 
sites  where  SAR  and  surface  roughness  measurements  were  obtained. 
In  addition,  the  ground  measurements  themselves  will  be  described. 

3.1  ERIM/CCRS  CV-580  SAR  SYSTEM 

From  1978  through  1984,  ERIM  operated  a  multi  frequency  and  mul¬ 
tipolarization  SAR  in  conjunction  with  CCRS.  This  SAR  system  was 
maintained  through  a  leasing  agreement  between  ERIM  and  CCRS  and  is 
called  the  CV-580  SAR  System.  This  facility  was  used  extensively 
over  the  past  six  years  in  a  variety  of  oceanographic  and  terrestrial 
research  programs.  The  leasing  agreement  between  ERIM  expired  at  the 
end  of  CY84.  At  this  time,  the  CV-580  SAR  System  is  not  available 
for  U.S.  Government  sponsored  research  programs.  Efforts  are  under¬ 
way  to  transfer  this  system  into  a  U.S. -owned  aircraft. 

The  CV-580  SAR  System  was  essentially  an  X-band  radar  to  which 
an  L-band  and  a  C-band  capability  was  added.  Because  of  this  con¬ 
figuration,  the  X-band  channels  were  always  available,  and,  in  ad¬ 
dition,  one  could  operate  either  the  L-band  or  the  C-band  channels 


simultaneously  with  the  X-band  channel.  For  each  frequency  band, 
two  orthogonal  polarizations  were  available.  The  radar  could  trans¬ 
mit  either  horizontal  or  vertical  polarization  and  receive  both  the 
parallel-  and  perpendicular-polarized  returns.  The  radar  antenna 
could  be  pointed  so  as  to  image  on  either  side  of  the  aircraft. 
Table  2  summarizes  the  relevant  radar  parameters  for  this  system. 

Figure  2  presents  a  functional  diagram  of  the  CV-580  SAR  System. 
The  major  components  of  the  SAR  are  its  antennas,  the  transmitters 
and  receivers,  the  calibration  signal  generator,  the  optical  film 
recorders,  the  digital  tape  recorder  and  the  X-band  real-time 
processor. 

Figure  3  presents  a  more  detailed  diagram  of  the  reference  oscil¬ 
lators,  transmitters  and  receivers  of  the  SAR  system.  This  diagram 
shows  how  the  L-band  and  C-band  wavelengths  are  converted  to  X-band. 

Figure  4  presents  a  diagram  of  the  calibration  signal  generator 
(CS6)  system,  which  is  used  as  a  monitoring  and  reference  source  for 
the  periodic  verification  of  the  SAR's  system  transfer  function. 
The  calibration  signal  generator  produces  synthetic  target  signals 
which  are  added  into  the  radar  receiver  (Walker  and  Larson,  1981; 
Larson,  et  al.,  1981;  1982;  1985).  The  form  of  the  synthetic  target 
corresponds  to  a  given  range,  and  the  intensity  is  controlled  by  an 
accurately  calibrated  R.F.  attenuator.  At  specific  times  during  the 
imaging  flight,  the  calibration  signals  are  inserted  into  the  re¬ 
ceiver  at  the  antenna  terminals,  so  that  the  radar  receiver  detects 
them  as  a  radar  return.  Using  a  sequence  of  such  signals  with  dif¬ 
ferent  intensities,  along  with  the  processed  signals  from  the  corner 
reflectors,  a  calibration  curve  can  be  produced  of  output  signal 
power  versus  radar  cross  section  (oq).  When  combined  with  measure¬ 
ments  of  precision  corner  reflectors,  CSG  measurements  can  be  used 
to  achieve  absolute  calibration  of  the  SAR  (see  Walker  and  Larson, 
1981;  or  Larson,  et  al.,  1985).  The  synthetic  targets  from  the  CSG 
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Figure  2.  Schematic  Diagram  of  ERIM/CCRS  X/C/L  Band  SAR  System 


' - - ►  L-BAND  CHANNEL 

Schematic  Diagram  of  CV-580  SAR  Receiver/Transmitter 


EXCITER 


AZIMUTH 

MOnULATlOH 


RANGE  X 
AZIMUTH 
GATE 


DIOITAl 

ATTruuATOR 


SYNTHETICl 

DOPPLER 


Synthetic 

APERTURE 
LENGTH,  AND 
RANGE  BANDUIOTH 


RANGE  DELAY 


CSG  PROGRAMMABLE 
MICROPROCESSOR 


AMpi iiunr 

CONTROL 


MODE  AUTO- PRt PROGRAMMED  SEQUENCE 

SELECTOR  MANUAL-USER  SUPPLIED  RANGE.  AMPLITUDE,  HAVUn.'TH 


SYNTHETIC 

TARGCTS 

TO  ptrriv! 

IliPljT', 


Figure  4.  Schematic  Diagram  of  CV-580  SAR  X/L  Calibration  Signal 
Generator 


L-band. 


Figure  5  presents  an  optically  processed  L-band  image  of  the 
calibration  signals.  Table  3  summarizes  the  signal  power  and  equiv¬ 
alent  radar  cross  sections  for  these  calibration  signals. 

The  three  SARs  which  comprised  the  CV-580  SAR  System  are  similar 
in  that  they  all  use  the  synthetic  aperture  technique  to  produce 
imagery  with  fine  cross-range  resolution.  These  SARs  used  pulse 
compression  to  achieve  fine  resolution  in  the  range  dimension.  Each 
recording  channel  was  adjusted  to  produce  imagery  of  a  selected  swath 
parallel  to  the  flight  direction.  The  width  of  the  imaged  swath  was 
determined  by  the  range  increment  of  signals  displayed  on  the  re¬ 
cording  system.  The  displacement  of  the  recorded  swath  from  the 
flight  line  was  adjustable  by  the  radar  operator.  The  CV-580  SAR 
System  is  extensively  described  by  Rawson,  et  al .  (1975)  and  CCRS 
(1983). 

The  SAR  phase  histories  were  recorded  onboard  the  aircraft  both 
optically  on  70  mm  film  and  digitally  on  high  density  digital  tape 
(HOOT).  The  optical  recording  system  recorded  four  data  channels, 
while  the  digital  system  was  restricted  to  recording  only  two  of  the 
four  data  channels  during  data  collection.  In  addition  to  the 
optical  and  digital  recorders,  a  real-time  digital  SAR  processor  on¬ 
board  the  CV-580  generated  X-band  imagery.  The  real-time  imagery 
was  not  intended  to  be  of  high  quality,  but  provided  "quick- look" 
data  necessary  for  system  performance  evaluation  and  planning  pur¬ 
poses  during  data  collection  missions.  It  was  also  used  to  select 
digital  data  of  interest  for  subsequent  processing  and  analysis. 
Table  4  presents  the  various  data  recording  options  available  on  the 
system.  Figure  6  schematically  illustrates  the  various  SAR  data  re¬ 
cording  and  processing  options. 
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X/L/C  Opticil  Processor 


Recording  and  Processing 


3.2  DATA  SETS 


There  are  three  sources  of  data  for  calibrated  L-band  terrain 
measurements  which  will  be  used  in  this  program:  (1)  data  collected 
during  a  US6S  program  in  April  1984;  (2)  data  collected  during  an 
ONR-sponsored  program  in  August/ September  1984;  and  (3)  L-band  SAR 
imagery  archived  at  ERIM.  These  data  sets  are  discussed  in  this 
section. 

3.2.1  USGS  SAR  DATA 

In  April  1984,  the  U.S.  Geological  Survey  sponsored  a  series  of 
SAR  data  collection  flights  near  Durham,  N.C.  In  conjunction  with 
these  flights,  a  series  of  test  passes  were  made  over  a  forested 
area  in  order  to  verify  the  operational  readiness  of  the  SAR.  This 
area  consisted  primarily  of  research  forests  owned  and  maintained  by 
Duke  University  and  is  collectively  known  as  the  Duke  Forest  (see 
Figure  7).  This  forest  consists  primarily  of  pure  and  mixed  1  to  50 
hectare  plots  of  pine  (southern.  Loblolly  and  Virginia),  whose  ages 
range  from  one  to  over  100  years  old.  There  are  also  mixed  stands 
of  hardwoods  (red  gum,  yellow  poplar,  river  birch,  sycamore,  beech 
and  white,  red,  black,  post  and  blackjack  oak)  scattered  throughout 
the  area. 

Duke  University  researchers  and  foresters  actively  manage  this 
forest  and  records  of  tree  species,  age  and  density  are  maintained 
for  each  plot  within  the  forest.  A  total  of  five  flights  were  made 
over  the  Duke  Forest  test  site  and  11  flights  over  other  sites  in 
this  area.  Both  horizontally-  and  vertically-polarized  imagery  were 
collected.  Table  5  summarizes  the  SAR  flights. 

During  both  the  USGS  and  ONR  flight  programs,  a  set  of  calibrated 
reference  reflectors  were  deployed  to  provide  calibration  references 
for  the  SAR  imagery.  The  types  of  reflectors  are  pictured  in  Figure 
8.  A  series  of  corner  reflectors  were  deployed  at  a  site  within 
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Duke  Forest  as  well  as  at  the  Asheboro  Airport.  The  position  of 
these  reflectors  are  presented  in  Figure  9,  while  their  sizes  are 
presented  in  Table  6. 


A  set  of  ground  measurements  were  obtained  from  14  separate  sites 
in  the  Duke  Forest.  These  sites  are  presented  in  Figures  10a  and 
10b  and  are  summarized  in  Table  7.  Data  collected  at  each  of  the 
selected  test  sites  consisted  of  surface  roughness  measurements, 
dielectric  constant  values,  soil  density  and  soil  moisture.  A  mea¬ 
sure  of  surface  roughness  was  obtained  by  "tracing"  the  surface  pro¬ 
file  on  graph  paper.  A  four  foot  section  of  chart  paper  was  taped 
to  a  plywood  board.  A  water-resistant  pen  was  then  attached  firmly 
to  a  wooden  rod  and  used  to  trace  the  surface  contour  on  the  chart 
paper.  Several  four-foot  long  tracings  of  the  surface  contour  were 
obtained  from  each  field.  A  photograph  illustrating  the  surface 
contour  tracing  method  is  presented  in  Figure  11. 

Soil  samples  were  obtained  from  each  test  site  and  used  to  deter¬ 
mine  soil  density,  moisture  content  and  dielectric  constant.  In 
addition,  a  portable  instrument  was  used  to  obtain  a  measurement  of 
dielectric  constant  on  site.  This  instrument  consists  of  a  "Q"  meter 
operating  at  100  MHz  and  a  coaxial  cavity  resonator  for  measurements 
at  L-band,  1200  fflz. 

3.2.2  ONR  DATA  COLLECTION 

In  August/September  1984,  the  Office  of  Naval  Research  sponsored 
a  series  of  data  collection  flights  near  and  over  Long  Island,  New 
York.  As  part  of  the  program,  a  series  of  24  passes  were  made  over 
a  calibration  array  located  a  Grumman's  Peconic  River  Airport  (see 
Figure  12)  during  eight  separate  missions.  The  terrain  surrounding 
the  Peconic  River  Airport  consists  primarily  of  agricultural  fields, 
but  also  contains  some  forested  land,  swamps,  and  coastal  marshes. 
Both  horizontally-  and  vertically-polarized  data  were  collected  over 
a  wide  variety  of  incidence  angles  (see  Table  8). 
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Asheboro  airport  calibration  array 
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Figure  9.  Reflector  Positions  at  Asheboro  and  Duke  Forest 
Calibration  Arrays 
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TABLE  6 

SUMMARY  OF  REFLECTORS  DEPLOYED 
DUKE  FOREST,  NORTH  CAROLINA 

AT  ASHBORO  AIRPORT  AND 
DURING  APRIL  1984 

Desiq 

Number  of  Reflector 

nation  Reflectors  Types 

Reflector 

Size 

(cm) 

X-Band 

(m^) 

L-Band 

(m^) 

A  ■  1  ■  "mi  —  ' 
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Site  Designation 
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TABLE  7 

DUKE  FOREST  TEST  SITES 


Site  Description 


Duke  Forest  (Blackwood  Division) 

Grass  Field 
Additional  Moisture 
10-15  -  16  percent 
15-20  -  15  percent 

Blackwood 

Rough  -  Short  Pine  -  Trees  -  Weeds 
Additional  Moisture  -  Different  Times 

No.  of  L-3  (Blackwood  Rd.) 

Cut  -  Rough  Area 

Fallen  Trees  -  One  Year  Growth 

LSG  Area 

Reflector  Site  -  Ashville  Airport 

Reflector  Site  -  L-3 

Forest  Sites  (Korstian  Division) 
30-Year  Lobolly  Pines 

Deciduous  Trees  -  30-50  Years 

26-50  Year  Lobolly  Pines 

70  Year  Lobolly  Pines 

20  Year  Lobolly  Pines 

Dense  Young  Pines 

Hardwoods  -  20  Feet  Tall 

20  Feet  -  Pines 

White-Red  Oak  -  60-80  Feet 


-  -p  .* 

%  ‘  /  -A*  -  1  ■ 


Location 


table  ft 

SUMMARY  OF  ONR  MULTICHANNEL  SAR  FI  lUriTS 


Frequency  an i 

•'u1  ar  1  /  it 

1  >n  - 

Digital 

Date 

Mission 

Pass 

* 

Area  ^ 

T  ime 

x-vy 

X-HH  x-'/H 

X-’IV  1 

-vv 

L-HH 

L-VH  t-HV 

26  August 

1 

1-7 

0 

1435-1544 

X 

X 

X 

X 

28  August 

2 

1-4 

0 

2044-2126 

X 

X 

X 

X 

30  August 

3 

1-2 

p 

1804-1828 

X 

X 

X 

X 

30  August 

3 

3-9 

0 

1837-2116 

X 

X 

X 

X 

3)  August 

4 

1-9 

0 

2026-2300 

X 

X 

X 

X 

31  August 

4 

1-12 

p 

2315-2357 

X 

X 

X 

X 

3  September 

5 

1-2 

p 

1726-1815 

X 

X 

X 

X 

5 

3-12 

0 

1832-21 W 

X 

X 

X 

X 

4  September 

6 

1-4 

p 

1853-2014 

X 

X 

X 

X 

6 

5-11 

0 

2029-2237 

X 

X 

X 

X 

5  September 

7 

1-3 

p 

1718-1802 

X 

X 

X 

X 

5  September 

7 

3-12 

0 

1821-2136 

X 

X 

X 

X 

6  September 

8 

1-2 

p 
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X 

X 

X 

X 

6  September 

8 

3-10 

0 

1304- 1614 

X 

X 

X 

X 
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9 

1-2 

p 
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X 

X 

X 

X 
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9 

3-6 

0 

1921-1951 

X 

X 

X 

X 

7  September 

10 

1-2 

p 

1405- 1502 

X 

X 

X 

X 

7  September 

10 

3-4 

p 

1517-1541 

X 

X 

X 

X 

7  September 

10 

5-8 

0 

1550-1657 

X 

X 

X 

X 

7  September 

II 

1- 

H 

2030-2036 

X 

X 

X 

X 

2- 

M 

2051-2055 

X 

X 

X 

X 

*  Areas 


P  -  Peconic  River  area 
0  -  Ocean  areas 
M  -  Mountain,  N.Y. 
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An  extensive  array  of  calibrated  reference  reflectors  were  de¬ 
ployed  at  the  Peconic  River  Airport.  The  positions  of  these  reflec¬ 
tors  are  presented  in  Figures  13  through  15,  and  the  size  of  these 
reflectors  are  summarized  in  Table  9. 

The  ground  measurements  discussed  in  the  previous  section  were 
also  collected  at  various  test  sites  out  and  around  the  Peconic  River 
Airport.  These  test  sites  are  listed  in  Table  10  and  are  identified 
on  the  diagrams  given  in  Figure  16. 

3.3  ARCHIVED  SAR  DATA 

There  is  an  extensive  set  of  calibrated  L-band  SAR  imagery  ar¬ 
chived  at  ERIM  which  is  available  for  use  in  this  program.  Of  par¬ 
ticular  interest  are  the  data  sets  listed  in  Table  11.  All  these 
data  sets  were  collected  after  1979,  when  the  CV-580  SAR  System 
underwent  an  extensive  system  modification  to  improve  the  system's 
calibration  capability.  All  the  imagery  from  the  missions  listed  in 
Table  11  contain  calibration  signals,  which  will  enable  calibration 
of  the  terrain  signatures  from  the  targets  of  interest. 

A  number  of  different  terrain  types  were  imaged  during  the  data 
collections  listed  in  Table  11.  These  experiments  are  discussed  in 
more  detail  by  de  Ritter,  et  al.  (1984),  Burns,  et  al .  (1984) 

Kasischke,  et  al.  (1983),  and  Shuchman,  et  al.  (1983).  Table  11 
summarizes  the  different  terrain  types  and  incidence  angles  which 
potentially  can  be  used  in  this  program  using  the  recently  collected 
and  archived  L-band  SAR  data. 


Figure  13.  Location  of  Calibrated  Reflector  Sites  at 
Peconic  River  Airport 
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Figure  14.  Location  of  Calibrated  Reflectors  at  Areas  1  and  2 
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TABLE  9 

SUMMARY  OF  CORNER  REFLECTORS  DEPLOYED  AT  PECONIC  RIVER 
AIRPORT  DURING  AUGUST/ SEPTEMBER  1984 


Designation 

Number  of 
Reflectors 

Reflector 

Type 

Reflector 

Size 

(cm) 

Radar  Cross 

Section 

X-Band 

(m^) 

L-Band 

(m2) 

A 

1 

AR1 

1538 

B 

3 

SQS 

119 

73828 

1429 

C 

1 

AR2 

— 

— 

631 

0 

2 

TRA 

119 

8200 

159 

E 

3 

TRA 

116 

7406 

143 

F 

2 

TRA 

91 

2805 

54 

G 

1 

TRA 

88 

2453 

48 

H 

2 

TRA 

68 

874 

17 

I 

3 

TRA 

61 

566 

11 

J 

1 

TRA 

52 

300 

6 

K 

1 

LLS 

— 

262 

— 

L 

2 

TRA 

45 

168 

3 

M 

3 

LLS 

— 

104 

— 

N 

1 

TRA 

39 

95 

2 

0 

1 

SQA 

22 

89 

2 

P 

2 

TRA 

38 

85 

1.5 

Q 

2 

LLS 

— 

32 

— 

R 

3 

TRA 

29 

30 

— 

S 

1 

SQA 

17 

30 

— 

T 

1 

TRA 

23 

11 

— 

U 

2 

TRA 

22 

10 

— 

V 

20 

SQS 

61 

5097 

99 

*AR1  -  Single  Polarization  Active  Radar  Calibration 
AR2  -  Dual  Polarization  Active  Radar  Calibration 
LLS  -  Luneberg  Lens 

SQA  -  Square  Aluminum  Precision  Corner  Reflector 

SQS  -  Square  Styrofoam  Corner  Reflector  (Aluminum  Coated) 

TRA  -  Triangular  Aluminum  Precision  Corner  Reflector 
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TABLE  10 

LONG  ISLAND  TEST  SITES 


Site  Description 


Strawberry  Field 
Orchard 
Smooth  Field 

Smooth  Field  -  Near  Sound  Rd. 

Potatoes  -  South  Side 

Potatoes  -  Row  N-S 
Row  N-S 

Row  Direction  -  345* 

Grass 

Corn  -  Stocks  -  Picked  but  Standing 
Dry  Stocks 

Cut  Corn 

Grass 

Cabbage 

Grass  Field  (VOR,  Reflectors) 

Grumman  Reflector  Site  -  4  Reflectors 

Grumman  Reflector  Site  -  8  Reflectors 

Grumman  Reflector  Site  -  26  Reflectors 

Short  Pine  Trees  -  Weeds 

(Similar  to  "West  of  L-3,  North  Carolina) 

Smooth  Field 

Smooth  Field 

Smooth  Field 

Smooth  Field 

Potatoes 

Smooth  Field 


Bare,  Smooth  Field 
N-S  Furrows 


Mountains 
Upper  N.T.  State 
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TABLE  11 

AVAILABLE  L-BAND  SAR  DATA  SETS 


Date  of 

Data  Collection 
June-July  1983,  84 
June-August  1983 

September  1980 
September  1980 


Location  -  Description 
_ of  Sites _ 

Sea  ice  and  ocean 

SE  Michigan,  urban,  rural 
areas,  Victoria,  B.C.; 
mountains,  ocean 

Michigan;  cultivated 
areas;  forests;  lakes 

Central  Canada  -  ground 
truth  available; 
cultivated  field  with 
caps,  base  fields,  tundra; 
forests;  deciduous  and 
coniferous  prairie  lands 


Incidence  Angle 
Range _ 

30‘-80* 

40* -85’ 


35-85 


40* -80* 


DATA  REDUCTION 


This  chapter  discusses  the  techniques  and  procedures  used  to  re¬ 
duce  the  SAR  and  ground  truth  data  collected  during  the  USGS  and  ONR 
experiments.  A  section  is  devoted  to  each  data  set. 

4.1  SAR  DATA 

The  SAR  data  recorded  during  the  two  missions  were  optically 
processed  into  imagery  using  ERIM's  tilted-plane  optical  processor 
(Kozma,  et  al.,  1972).  An  engineering  assessment  was  performed  on 
the  data  to  determine  which  passes  met  minimum  image  quality  stan¬ 
dards.  From  the  optically-processed  imagery,  specific  scenes  were 
selected  and  digitally-processed  into  imagery  using  ERIM's  Digital 
Strip  Map  Processing  facility  (Jackson,  et  al.,  1974). 

Figure  17a  and  17b  present  representative  examples  of  optically- 
processed  L-band  imagery  collected  during  the  USGS  flights.  These 
images  include  the  portion  of  the  Duke  Forest  where  the  corner  re¬ 
flectors  were  deployed.  Figures  18a  and  b  present  L-band  SAR  imagery 
of  the  Asheboro  Airport  Calibration  Array. 

Figures  19a,  b  and  c  present  digitally-processed  L-band  imagery 
of  the  Peconic  River  Airport  gathered  during  the  ONR  data  collection 
mission.  These  images  contain  the  corner  reflector  array  as  well  as 
the  ground  test  sites  (see  Figures  13  and  16).  Figures  20a  and  20b 
present  optically-processed  L-band  imagery  collected  over  a  forested 
mountain  region  in  upper  New  York  for  horizontal  and  vertical  polari¬ 
zation,  respectively. 

In  the  remainder  of  this  section,  we  will  discuss  the  techniques 
used  to  calibrate  the  digital  L-band  SAR  data  and  present  examples 
of  calibrated  terrain  measurements. 


Figure  17b.  Optically-Processed  L-Band  (VV)  SAR  Imagery  of 
Duke  Forest  Test  Site  (USGS-l,  Pass  l,  April  8 
198A;  Arrow  Indicates  Calibration  Array  Locati 


Figure  1?<b.  Upt ically-Processeo  L-Hand  (VV)  SAR  Im, 
Asheroro  Airport  Calibration  Array  (US' 
April  x,  1984;  Arrows  Indicate  Array  Ij 


Figure  19a.  Digitally-Processed  L-Band  (HH)  SAR  Imagery  of 
Peconic  River  Calibration  Array  (SARSEX-10, 
Pass  1,  September  7,  1989;  Arrows  Indicate 
Reflector  Location) 


Figure  19b.  Digitally-Processed  L-Band  (VV)  SAR  Imagery  of 
Peconic  River  Calibration  Array  (SARSEX-10, 
Pass  9,  September  7,  1989;  Arrows  Indicate 
Reflector  Location) 


Figupe  19c.  Digitally-Processed  L-Band  (HH)  SAR  Imagery  of 

Ground  Measurement  Sites  Near  Peconic  River  Airport 
(SARSEX-10,  Pass  2,  September  7,  1989) 


Figure  20a.  Optically-Processed  L-Band  (HH)  SAR  Imagery  of  Upper 
New  York  (SARSEX-IL  Pass  1,  September  7,  1984 


4- 


E 

5 


4.1.1  CALIBRATION  OF  SAR  DATA 


The  procedure  to  convert  the  recorded  SAR  image  intensities  into 
radar  cross-section  actually  involves  two  steps:  image  calibration 
and  radiometric  corrections.  This  procedure  is  carried  out  in  the 
slant  range  plane.  If  a  calibrated  image  in  ground  plane  coordinates 
is  desired,  a  slant-to-ground  range  conversion  can  be  subsequently 
performed.  The  algorithms  to  accomplish  calibration  are  described 
in  the  following  paragraphs. 

I ma g^e  C a  1  i  b r at  ion 

The  image  calibration  algorithm  translates  the  image  intensity 
values  into  actual  radar  cross-sections.  This  is  accomplished  by 
comparing  the  image  intensities  of  the  SAR  data  to  be  calibrated  with 
known  image  intensities  produced  by  the  "reference"  calibrated 
reflectors. 

The  calibration  of  a  SAR  image  begins  with  the  well  known  radar 
return  equation 

HP,G^(e)xo 

”  ■  ---TT- *  %  l‘> 

4it  R 

where  P  is  the  measured  image  power, 
is  the  system  noise  power, 

R  is  the  slant  range, 

X  is  SAR  wavelength, 

P^  is  the  transmitter  power, 

G(e)  is  the  antenna  gain, 
e  is  the  incidence  angle, 
o  is  the  radar  cross-section,  and 
H  is  the  SAR  gain  from  receiver  input  to  image. 


3  4 

An  R  term  is  used  rather  than  an  R  term  because  a  1/R  fac¬ 
tor  is  applied  to  the  image  during  the  azimuth  compression  stage  of 
SAR  image  processing.  Equation  (1)  can  be  rewritten  as 


(P  -  P^)(4n)¥ 
HxP^G^(g) 


(2) 


If  a  SAR  images  a  scene  with  a  series  of  targets  with  known  radar 
cross-sections  of  different  values,  then  a  relationship  between  o 
and  P  can  be  established,  and  can  be  calculated,  provided  the 
transmitted  power  (P^)  and  antenna  gain  (G(e))  are  known 

The  procedure  used  to  calibrate  aircraft  SAR  imagery  is  to  relate 
the  data  to  be  calibrated  (o.,  where  the  i  subscript  denotes  the 
scene  to  be  calibrated)  to  the  reference  scene  radar  cross-sections 
(o^,  where  the  r  subscript  denotes  the  reference  scene),  i.e.. 


-  ’’hO 

V  ■  (P^  - 


or 


-  '’ni>  - 


: _ vKlH 


(4) 


The  calibration  procedure  determines  all  the  terms  in  Eq.  (4), 

and  thus  allows  for  calculation  of  for  any  P^  in  the  scene. 

In  Eq.  (4),  all  the  terms  are  known  or  can  be  calculated  except  for 

0  IIP  -  P.,  )  and  The  relationship  between 

r  '  r  Nr  i  r  r 

and  P^  is  determined  by  measuring  the  intensities  of  a  set  of 
calibrated  retroref lectors  deployed  in  the  calibration  array  and 


52 


imaged  by  the  SAR.  The  measured  power  intensity  (P^)  values  are 


then  regressed  against  the  radar  cross-section  (a^)  values, 


Pp  =  a  0^  +  b 


(5) 


where  a  is  the  linear  regression  coefficient  and  b  is  the 
y-intercept.  If  we  assume  the  y-intercept,  b,  is  equivalent  to  the 


system  noise,  then  Eq.  (5)  can  be  rewritten  as 


a  = 


-  V) 


(6) 


or. 


1 

a  =  (P’-TTT  ^  ’ 

'  r  Nr^ 

The  next  step  in  the  calibration  procedure  is  to  calculate 
and  H..  Using  Eqs.  (1)  and  (6),  can  be  calculated.  H^.  is 
determined  by  using  the  intensities  of  the  synthetic  targets  gener¬ 
ated  by  the  calibration  signal  generator.  The  CS6  produces  a  series 
of  synthetic  point  targets,  each  target  with  a  different  attenuation 
setting.  The  CSG  signal  is  the  same  transmitted  pulse  as  used  for 
the  data.  Therefore,  the  CSG  signals  also  monitor  the  transmitted 
waveform.  The  calibration  signals  are  inserted  at  the  antenna  prior 
to  and  after  each  data  collection  pass,  and  are  recorded  on  the  dig¬ 
ital  tape  in  the  same  fashion  as  the  SAR  phase  histories.  They  are 

then  processed  into  images.  Thus  any  differences  in  the  intensity 
in  the  processed  image  of  a  specific  signal  between  passes  represents 
a  change  in  the  SAR  system  gain  (H). 

The  power  intensity  of  the  calibration  signals  generated  after 
the  calibration  reference  pass  and  the  image  to  be  calibrated  are 
measured.  These  intensities  are  then  plotted  on  a  curve,  using  the 

input  power  from  the  CSG  as  the  axis.  This  plotting  results  in  two 

straight  lines  (one  for  each  pass),  each  with  the  same  slope.  The 


vertical,  relative  difference  between  these  two  lines  represents  the 
amount  the  system  gain  for  the  reference  pass  (H^)  has  to  be  ad¬ 
justed  to  obtain  a  system  gain  for  the  image  to  be  calibrated  (H.). 

In  the  implementation  of  Eq.  (4)  to  calculate  radar  cross- 
sections,  logarithms  of  the  values  are  usually  taken,  and  all  plots 
and  computations  are  carried  out  in  dB  for  convenient  numerical 
evaluation. 

4.1.2  INITIAL  RESULTS 

To  date,  two  calibration  data  sets  have  been  processed  and  cali¬ 
bration  relationships  obtained.  Examples  of  the  calibration  rela¬ 
tionships  are  given  in  Figure  21.  Figure  21  presents  the  L-band 
calibration  relationship  obtained  for  the  Long  Island  flights  of  31 
August  and  7  September  1984  using  the  calibrated  reflectors  to  pro¬ 
vide  the  reference  for  calibration  RCS  values. 

As  initial  verification  of  the  L-band  calibration  algorithm,  the 
SAR  data  obtained  on  7  September  1984  was  used  to  calibrate  the  SAR 
data  of  31  August  1984.  Both  SAR  data  sets  include  the  reflector 
array  at  the  Peconic  Airport,  Long  Island.  The  values  for  each  re¬ 
flector  of  the  31  August  data  were  determined  using  the  7  September 
calibration  curve.  These  values  were  then  compared  to  the  known 
nominal  values  for  each  reference  reflector;  the  results  are  pre¬ 
sented  in  Figure  22.  A  standard  deviation  of  about  3  dB  is  obtained 
using  values  from  21  reference  reflectors. 

Using  the  calibration  relationship,  the  values  of  for  sev¬ 
eral  of  the  fields  selected  for  initial  analysis  were  determined. 
Values  of  Oq  are  given  in  Table  12  for  several  test  sites.  These 
values  are  obtained  by  evaluation  of  a  histogram  of  the  radar  cross 
sections  for  these  sites.  The  histograms  of  RCS  values  for  the  area 
near  the  reference  reflector  array  is  given  in  Figure  23;  other  RCS 
histograms  are  included  in  the  Appendix. 
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Figure  22.  Plot  of  Predicted  versus  Actual  Radar  Cross 
Section  Values  for  L-Band  Comer  Reflectors 


TABLE  12 

VALUES  OF  RADAR  CROSS  SECTION  (on)  FOR 

SEVERAL  TEST 

SITES  IN  THE 

PECONIC 

RIVER  AIRPORT 

TEST  AREA 

(L-Band  Horizontal , 

Polarization, 

63.9*  Incidence  Angle) 

F  ield 

a 

7  * 

(dBm^) 

% 

0 

0 

(S.D 

Forest  Area 

Near  Reflector  Site 

1 

-16.6 

0.012 

O.C 

2 

-15.0 

0.014 

o.c 

3 

-14.2 

0.017 

O.C 

Grass  RF2 

-23.0 

0.002 

O.C 

Grass  RF3 

-19.9 

0.004 

o.c 

Grass  RF3 

-20.4 

0.004 

o.c 

E2 

-12.3 

0.023 

o.c 

El 

-11.0 

0.029 

o.c 

FI 

F2 

-17.6 

0.008 

o.c 

FI 

0.008 

o.c 

G2 

-14.9 

0.015 

o.c 

G1 

-11.9 

0.003 

o.c 

*Total  Radar  Cross  Section 

is  in  dB  with 

reference  to 

one  square 

meter  for  one  pixel.  Pixel  size  = 

4.2  GROUND  DATA 


The  soil  samples  collected  during  the  SAR  overflights  were  used 
to  obtain  three  parameters:  dielectric  constant,  soil  moisture  and 
soil  density.  Dielectric  constant  measurements  were  made  using  a 
coax  cavity  at  L-band.  Estimates  of  the  precision  of  the  dielectric 
measurements  were  obtained  by  measuring  different  sub-samples  of  the 
same  soil  sample.  The  standard  deviation  obtained  was  about  5  per¬ 
cent  of  the  mean.  Soil  moisture  measurements  were  obtained  in  the 
laboratory  using  samples  collected  at  each  site.  A  known  volume  of 
each  sample  was  weighed,  then  heated  at  300"F  in  an  oven  for  a  period 
of  at  least  8  hours.  The  dryed  sample  was  then  weighed.  The  before 
and  after  weights  were  used  to  compute  moisture  content  by  weight 
and  by  volume.  Soil  density  measurements  were  obtained  by  weighing 
a  known  volume  of  each  soil  sample.  The  measured  soil  parameters 
are  summarized  in  Table  13  for  the  Duke  Forest  test  sites  and  in 
Table  14  for  the  Peconic  River  Airport  test  sites. 

The  raw  surface  profile  measurements  were  digitized.  Figure  24 
presents  an  example  of  a  digitized  height  profile  obtained  from  the 
raw  measurements.  This  profile  was  obtained  from  the  grass  field 
near  reflector  Site  3  at  the  Peconic  River  Airport.  These  data  will 
be  used  as  an  example  throughtout  the  remainder  of  this  section. 
Examples  of  these  raw  measurements  for  each  of  the  fields  is  included 
in  the  Appendix.  These  raw  measurements  were  then  used  to  calculate 
a  set  of  surface  descriptors. 

The  surface  height  profiles  to  be  correlated  with  the  observed 
scattering  patterns  are  considered  as  statistical  samples  of  a  random 
process.  Surface  height,  h(x),  is  considered  as  a  random  variable 
and  various  statistical  measures  are  derived.  These  measures  will 
be  used  in  attempts  to  correlate  surface  properties  with  observed 
scattering  patterns.  At  this  time  only  horizontally-polarized  data 
over  a  very  limited  range  of  incident  angles  have  been  considered. 
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3  (Reflector  Array)  5-10  cm  1.04  1^.93  1.05  3.69  0.?90 

3  (Grassy  Area)  10-12  cm  0.99  8.39  0.99  3.40  0.125 


Therefore,  initial  comparisions  are  made  between  the  statistics  of 
the  SAR  reflectivety  data  for  several  fields  but  at  one  incident 
angle.  The  next  year's  effort  will  expand  the  analysis  to  obtain 
correlation  if  possible,  over  a  range  of  incident  angles  and  for  both 
H  and  V  polarization. 

Using  the  surface  height  profiles,  additional  roughness  proper¬ 
ties  were  calculated: 

1.  Spectrum  of  the  surface  height; 

2.  Slope  Distribution; 

3.  Radius  of  Curvature;  and 

4.  Auto-Correlation  of  the  surface  heights. 

Surface  height  spectra  are  obtained  by  taking  an  FFT  of  the  sur¬ 
face  height  profiles.  An  example  of  the  spectra  obtained  for  the 
height  profile  in  Figure  24  is  given  as  Figure  25.  Additional  spec¬ 
tra  are  given  in  the  Appendix. 

Surface  slope  information  Is  obtained  by  differentiation  of  the 
surface  height  profiles.  The  histogram  for  the  calculated  slope 
values  from  the  example  field  is  given  as  Figure  26.  Values  of  the 
mean  slope  for  a  number  of  the  test  areas  are  listed  in  Table  15. 

Values  for  slope  spectra  S^(k)  can  also  be  obtained  directly 
from  the  FFT  of  the  slope  function  or  from  the  spectra  of  the  height 
profile  data  using  the  relationship: 

SjCk)  =  K\(k)  (8) 

where  S^(k)  =  surface  height  profile  spectrum 
S^(k)  =  slope  spectrum 

An  example  of  a  surface  slope  spectra  is  given  in  Figure  27. 

An  additional  parameter  that  is  related  to  the  scattering  cross- 
section  of  surfaces  is  the  radii  of  curvature  of  surface  structure. 
The  initial  approach  at  obtaining  radii  of  curvature  statistics  from 
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SURFACK  HEEGHT  SPECTRA  —  LONG  ISLAND 


Figure  25.  Example  of  the  FFT  Spectra  of  the  Surface  Height 
Profile 


Figure  26.  Example  of  the  Histogram  of  Surface  Slope  Values 


Surface  Slope  Spectra 


the  measured  surface  height  profile  is  to  utilize  the  generated  curve 
of  slope  vs.  distance  and  calculate  second  derivatives.  Thus,  for  each 
sample  point  a  radius  of  curvature  is  calculated  using  the  stan¬ 
dard  relationship 


.  (1  * 

*^c  h” 


(9) 


where  h'  =  surface  slope  [1st  derivative  of  h(x)] 
h"  =  second  derivative 

Examples  of  mean  value  of  radii  of  curvature  obtained  are  given  in 
Table  16  and  of  the  histogram  for  the  radius  of  curvature  in 
Figure  28. 

The  autocorrelation  function  of  the  surface  profile  is  defined  as 


h(x’)h(x  -  x^)  dx  (10) 

h^  - 

A  basic  computer  program  performs  the  correlation  with  an  example 
autocorrelation  function  given  in  Figure  29  and  in  the  Appendix. 

Additional  histograms  for  both  slope  distributions  and  radii  of 
curvature  were  also  generated  and  examples  of  each  are  given  in  the 
Appendix. 


TABLE  16 

SCATTERING  COEFFICIENTS  OBTAINED  FROM  REFLECTOR  SITE  3 


Date 

31  August 
7  September 


Incident  Angle 
(degrees) 

56.5 

63.9 


0.024 

0.004 


Figure  28.  Example  of  Histogram  of  Radii  of  Curvature  Values 


139.75  209.13  278.5  347.88  417.25  486.63  556.00 

Sample  points  —  divide  by  4  to  obtain  distance  (inches) 

Example  of  An  Auto  Correlation  Surface  Roughness 
Profile 


5 

DATA  ANALYSIS 


Descriptions  of  backscatter  from  terrain  and  ocean  are  required 
by  both  the  microwave  remote  sensing  researcher,  particularly  those 
utilizing  SAR  imagery  and  by  the  radar  system  design  engineers.  It 
is  interesting  to  point  out  that  for  the  remote  sensor  scientist 
backscatter  forms  the  desired  signal.  And  by  contrast,  for  the  en¬ 
gineers  designing  detection  radars  surface  scattering  or  clutter  be¬ 
comes  the  noise  competing  with  desired  signals.  The  important  radar 
parameters  for  which  backscatter  characteristics  will  vary  include 
wavelength,  aspect  angle,  polarization  and  resolution.  There  is  a 
very  large  set  of  terrain  types  and  conditions  that  are  encountered. 
These  range  from  rough  to  smooth  surfaces,  wet  to  dry  conditions  and 
bare  surfaces  and  vegetated  areas. 

Descriptions  of  the  radar  sensor  parameters  can  be  quantified; 
however,  quantifying  the  terrain  scattering  surfaces  has  proven  to 
be  very  difficult.  Several  measures  have  been  used  to  describe  ter¬ 
rain.  For  example,  a  measure  of  surface  roughness,  dielectric  con¬ 
stant,  density  and  moisture  content,  height  and  density  of  vegetation 
and  others.  Qualitative  descriptions  such  as  photographs,  soil  type, 
vegetation  type  are  also  used  to  classify  terrains. 

During  the  past  30  years  or  so,  various  researchers  have  studied 
the  scattering  characteristics  of  terrain.  An  extensive  research 
effort  conducted  by  Peake  (1960)  during  the  late  1950's  to  1965  pro¬ 
vided  the  first  systematic  approach  to  the  study  of  scattering  prop¬ 
erties.  This  work  consisted  of  analysis  and  measurement  of  the  angu¬ 
lar  response  of  scattering  from  various  terrain  types  and  surfaces. 
Also,  the  dependency  on  wavelength  was  considered  for  both  horizontal 
and  vertical  polarization.  Subsequent  work  from  the  mid  1960's  ex¬ 
panded  on  the  analysis  with  many  publications  appearing.  An  excel¬ 
lent  summary  is  given  in  Beckmann  and  Spizzichino  (1963).  A  more 
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recent  compilation  of  scattering  theories  is  given  in  a  two  volume 
series  by  Ishimaru  (1978).  Also,  an  excellent  summary  giving  results 
from  a  number  of  scattering  measurement  programs  is  given  by  Long 
(1983).  And  a  summary  of  recent  work  on  the  electromagnetic  interac¬ 
tion  problem  is  given  by  Ulaby,  et  al.  (1981  and  1982). 

All  of  this  effort  on  electromagnetic  scattering  has,  in  one  way 
or  another,  attempted  to  relate  a  scattering  model  to  parameters 
describing  various  surfaces.  Radar  systems  engineers  have  developed 
models  to  describe  the  radar  clutter  environment.  Clutter  models 
have  been  developed  using  results  obtained  from  the  various  measure¬ 
ments  and  analysis  described.  Several  researchers  have  reported  on 
the  development  of  clutter  models  such  as,  for  example  Greenstein, 
et  al.  (1969). 

The  availability  of  data  from  calibrated  synthetic  aperture  radar 
(SAR)  systems  has  provided  an  additional  source  of  scattering  data. 
Also,  the  SAR  has  been  used  in  remote  sensing  collection  demonstra¬ 
tion  programs  to  obtain  data  from  very  large  areal  areas  and  from 
remote  locations  (MIZEX  1983-84).  The  calibrated  SAR  has  also  been 
utilized  in  the  study  of  clutter  characteristics  (Larson,  et  al., 
1981). 

5.1  SCATTERING  MODELS 

An  objective  of  this  program  is  to  utilize  the  measured  clutter 
data  in  scattering  models  in  order  to  improve  the  capability  of  pre¬ 
dicting  scattering  properties.  Analysis  during  the  first  year  has 
been  directed  toward  the  selection  of  models  to  describe  scattering 
in  terms  of  the  surface  properties  of  the  various  test  fields  using 
an  example  of  a  statistical  method  of  describing  the  measured  distri¬ 
butions  of  the  backscatter  coefficients.  For  the  initial  analysis, 
data  has  been  selected  from  several  test  sites  in  the  Long  Island 
area.  The  results  obtained  from  analysis  during  the  first  year  are 
described  in  this  section. 


Results  from  analysis  of  the  surface  measurement  data  indicate 
that  most  of  the  test  sites  satisfy  the  slightly  rough  or  rough  sur¬ 
face  criteria  for  23  cm  wavelength  data.  A  number  of  semi-empirical 
models  for  the  scattering  cross-section  from  rough  surfaces  have  been 
proposed  (Ruck  1970).  Several  of  these  models  include  incidence 
angle,  wavelength  and  some  measure  of  surface  roughness  properties 
as  input  parameters.  One  of  the  models  listed  by  Ruck,  proposed  by 
Spetner  and  Katz  (1960),  has  been  selected  as  a  candidate  for  use  in 
relating  surface  properties  to  the  measured  scattering  pattern.  The 
model  is  given  as 


^-tan2,,/2s2 


(11a) 


(Xq  >  x^)  (11b) 


where,  K^,  Kg,  Kg  are  constants  defined  in  the  reference  in 
terms  of  the  surface  roughness  properties 

Xq  -  wavelength 

x^  -  upper  limit  on  slope  spectrum 
2 

Sq  -  mean-square  slope  of  surface 

Barrick  and  Peake  (1967)  have  derived  a  scattering  model  for  a 
slightly  rough  surface  given  as 

a  =  4i?  h^  cos"^  ^  W(-2k  sin  o.,  0)  (12) 

pq  0  1  pq  1  '  '  ' 

where,  for  backscatter 


kjj  -  wave  number 


2 

h  -  rms  surface  height 

Opq  -  scattering  matrix  elements  for  p,  q  polarization 
W  -  surface  height  power  spectrum 

(13) 

For  situations  where  very  rough  surface  conditions  are  applica¬ 
ble,  the  Barrick  (1968)  specular  point  model  will  be  applied.  The 
form  of  this  model  is  given  as 

o  =  IT  sec%  p(s)R(o)^  (14) 

where  y  -  local  angle  of  incidence 

p(s)  -  slope  probability  density  function 
R(o)  -  Fresnel  coefficient,  normal  incidence 

Selection  of  the  candidate  models  described  is  tentative  and  is 
based  on  the  use  of  surface  descriptions  that  can  be  obtained  from 
emperical  data.  It  is  anticipated  that,  during  the  next  year  effort 
models  will  be  modified  to  best  fit  the  data  available.  Also,  other 
models  will  be  investigated  for  applicability  such  as  that  developed 
recently  by  Winebrenner  and  Ishimaru  (1985). 

Although,  thus  far  only  data  obtained  at  one  incidence  angle  have 
been  used  in  the  analysis,  a  comparison  of  values  of  different 
angles  from  one  site  near  the  reflector  array  has  been  made.  This 
site  is  a  grass  area  near  reflector  site  3  that  is  included  on  both 
the  31  August  and  7  September  data.  The  values  of  obtained  are 
given  in  Table  16,  both  for  horizontal  polarization.  Two  val¬ 
ues  obtained  7  days  apart  for  the  same  location  are  seen  to  differ 
by  7.4  dB.  Field  roughness  data  indicated  that  the  field  was  slight¬ 
ly  rough  for  the;  L-band  wavelength  of  23  cm.  Using  the  Barrick 
(1967)  model  for  slightly  rough  surfaces,  given  in  Eq.  (12),  a  ratio 
of  Oq  values  for  the  two  measurements  was  obtained. 


Ground  truth  records  showed  that  the  area  was  very  wet  during 
the  time  of  data  collection  on  31  August,  due  to  heavy  rain  the  pre- 
ceeding  morning.  Although  dielectric  constant  measurements  were  not 
made  on  31  August,  the  fact  that  the  field  was  saturated  would  put 
the  real  part  of  the  dielectric  constant  over  10  with  a  very  large 
loss  tangent.  Using  these  estimated  values  and  the  measured  dielec¬ 
tric  constant  values  obtained  on  7  September  to  calculate  values  for 
“hh*  ®  value  for  the  a  ratio  of  about  7  dB  is  obtained.  Note  that 
there  is  a  7.4  degree  difference  in  incidence  angle  for  the  two  data 
sets. 

Additional  calibration  data  points  will  be  obtained  during  the 
next  year's  effort  to  refine  the  calibration  and  to  extend  the  cali¬ 
bration  to  include  the  North  Carolina  data  and  to  the  Mountain  data 
of  7  September. 

5.2  SAR  DATA  CHARACTERIZATION 

Since  the  attributes  which  characterize  most  of  the  earth's 
natural  surfaces  (such  as  terrain,  water,  or  ice)  are  essentially 
random  in  nature,  it  is  necessary  to  describe  the  scattering  from 
such  surfaces  in  terms  of  statistical  parameters.  The  ideas  to  be 
suggested  and  examined  in  this  section  proceed  somewhat  in  the  re¬ 
verse  of  the  usual  treatments.  That  is,  we  hypothesize  that  there 
is  at  least  one  probability  distribution  in  terms  of  which  surface 
scattering  can  be  well  described.  We  further  expect  that  the  param¬ 
eters  of  such  a  distribution  can  be  used  to  separate  not  only  types 
of  surfaces  but  also  incidence  angles  (and  perhaps  polarizations) 
for  which  a  particular  surface  is  examined. 

The  distribution  suggested  fcr  this  role  is  the  inverse  Gaussian 
(16)  distribution.  This  distribution  has  been  around  since  the  days 
of  Schroedinger  (1915)  who  used  it  in  the  study  of  first  passage  time 
in  Brownian  motion.  Tweedie  (1957)  proposed  the  inverse  Gcussian 


name  because  of  the  inverse  relationship  which  exists  between  the 
cumulant  generating  function  of  the  first  passage  time  distribution 
and  that  of  the  Gaussian  distribution.  Excellent  discussions  of  this 
distribution  in  terms  of  its  history,  its  properties,  and  some  of 
its  applications  are  given  by  Folks  and  Chhikara  (1978)  and  by  Cheng 
and  Amin  ( 1981 ) . 

Tweedie  suggests  at  least  four  forms  for  the  IG  density  function. 
The  one  which  we  shall  adopt  for  our  purposes  is  given  by 


f  (x;  a,  li,  \)  = 


2jr(x  -  a) 


x(x  -  g  -  ii)‘ 

2ti^{x  -  a) 


where  u  is  the  mean  of  the  population,  x  is  equal  to  the  square-root 

1/2 

of  one-half  the  mean,  i.e.,  (u/2)  ,  and  x  is  defined  as 

3 

x=^  (16) 

S 

where  S  is  the  standard  deviation  of  population. 

One  observes  immediately  that  the  density  function  depends  upon 
the  three  parameters:  a,  u,  and  x.  It  is  proposed  that,  in  the 
space  of  these  three  parameters,  points  corresponding  to  samples  from 
different  terrain  types  will  be  well  separated.  And  not  only  will 
points  for  different  terrain  types  be  quite  distinct,  but  also  dif¬ 
ferent  regions  of  the  three-parameter  space  will  characterize  differ¬ 
ent  surface  types.  The  reasons  for  believing  that  these  ideas  hold 
some  validity  is  that  they  have  been  partially  tested  on  SAR  ice  data 
with  promising  results  (Maffett,  1979).  Eventually  we  wish  to  test 
the  fit  of  the  IG  distribution  to  the  cumulative  histograms  for  vari¬ 
ous  terrain  samples  by  some  standard  nonparametric  test  such  as,  for 
example,  a  Komologoroff-Smirnov  test  (see,  for  example.  Folks  and 
Chhikara,  1978  or  Conover,  1980). 


Our  reasons  for  choosing  the  approach  outlined  above  and  for 
choosing,  in  particular,  the  16  distribution  are  based  on  the  long¬ 
standing  and  rather  successful  application  of  the  log-normal  distri¬ 
bution  to  describe  clutter  data  of  various  sorts.  It  is  planned  to 
compare  calculated  IG  fits  with  log-normal  fits  for  the  same  data 
samples.  We  believe  that  the  16  distribution  will  offer  an  improve¬ 
ment  over  the  log-normal  because  of  its  tail  behavior.  The  tail  of 
the  IG  distribution  goes  more  quickly  to  zero  than  does  the  tail  of 
the  log-normal  distribution.  This  fact  should  be  advantageous  in 
describing  terrain  (or  other  earth  surface)  data  since  there  are  no 
infinite  values  of  clutter  RCS.  There  exists  exact  sampling  distri¬ 
bution  theory  for  the  IG  distribution  and  it  is  a  reasonably  easy 
distribution  to  work  with.  These  facts  make  the  IG  distribution  an 
attractive  candidate  to  describe  skewed  terrain  data. 

SAR  images  used  in  this  research  are  recorded  and  processed  to 
give  data  which  are  proportional  to  square  root  of  relative  power. 
Since  it  may  be  desirable  to  examine  resolution  cells  of  different 
sizes,  each  datum  of  a  sample  is  squared  before  any  analysis  is  car¬ 
ried  out.  In  this  way,  resolution  cell  size  may  be  varied  and  the 
results  normalized  to  the  resolution  cell  area.  The  result  is  nor¬ 
malized  radar  cross-section  (NRCS),  "o."  One  problem  to  be  investi¬ 
gated  is  the  effect  of  resolution  cell  size  on  the  distribution  from 
which  the  sample  has  been  drawn. 

The  chosen  sample  is  now  used  to  find  estimates  (a,  u,  and  x) 
for  the  parameters  a,  u,  and  x  of  an  IG  distribution.  The  estimates 
used  are  maximum  liklihood  (ML)  ones  given  by  Cheng  and  Amin  (1981) 
and  are  numerically  straightforward.  Their  formulas  require  an  iter¬ 
ation,  but  in  actual  practice  (and  especially  for  small  samples)  it 
is  almost  never  necessary  to  go  beyond  a^.  The  estimates  a,  u, 
and  X  are  calculated  as: 
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where  x^  is  the  smallest  datum  of  the  sample  and  x  and  s  are  the 
sample  mean  and  standard  deviation,  respectively. 

To  illustrate  the  nature  of  these  parameters,  an  example  will  be 
presented.  A  small  sample  (n  =  256  points)  from  a  SAR  image  of  a 
grassy  field  was  selected.  Each  datum  is  squared  and  considered  in¬ 
dividually;  that  is  to  say,  each  squared  datum  is  proportional  to 
the  power  returned  from  a  resolution  cell  of  dimensions  1.5  m  by 
1.44  m.  In  this  example  only,  the  datum  is  not  divided  by  the  reso¬ 
lution  cell  area.  The  smallest  squared  datum  had  a  value  4,  the 
largest  had  a  value  37636.  The  mean  was  5639.918  and  the  standard 
deviation  5651.794.  Using  the  above  formulas,  we  fird  the  ML  esti¬ 
mates  for  a,  u,  and  x  to  be  -1159.568,  6799.486,  and  8658.059,  re¬ 
spectively.  Only  the  first  step  of  the  iteration  was  used  since  the 
second  step  a  was  identical  with  the  first.  Figure  30  presents  a 
histogram  of  the  data  and  the  IG  density  resulting  from  the  above 
parameter  estimates.  The  fit  seems  quite  reasonable  considering  the 
small  sample  size. 


Figure  30.  Characterization  of  SAR  Data  using  the  Inverse 
Gaussian  Distribution 


5.3  EXAMPLE 


A  preliminary  examination  of  the  IG  distribution  inputs  (a,  u, 
and  x)  from  12  extracts  from  SAR  images  has  been  started.  These  ex¬ 
tracts  are  comprised  of  samples  from  3  grass-clutter  areas,  6  field 
areas,  and  3  forest  areas.  The  digital  output  from  each  of  these 
areas  has  been  described  above.  In  the  processing  to  be  discussed 
here,  the  data  have  not  been  converted  to  by  dividing  sums  of 

squares  pixel  outputs  by  appropriate  areas.  If  various  sampling 

methods  are  to  be  compared,  this  normalization  should  be  carried  out. 

Samples  of  sizes  ranging  from  900  to  12800  pixels  have  been 

processed  so  as  to  vary  the  effective  resolution  by  grouping  various 
sized  sets  of  pixels  (for  example,  1  x  1,  3  x  3,  and  so  on).  Esti¬ 
mates  for  the  IG  parameters  a,  u,  and  x  have  been  computed  for  each 
of  the  various  samples.  Table  17  summarizes  these  results  as  well 
as  the  means  and  standard  deviations  of  the  samples.  To  some  degree, 
the  'X,  u,  and  x  parameter  values  duplicate  some  of  the  information 
supplied  by  the  mean  and  standard  deviation.  For  example,  consider 

Field  FI  compared  with  Forest  3.  The  means  and  standard  deviations, 
as  well  as  the  u  values  compare  favorably.  The  a  and  x  values,  how¬ 
ever,  even  though  similar,  are  distinct  enough  to  separate  these  two 
scenes,  y  and  x  are  only  partially  interpretable  as  location  and 
scale  parameters  since  (if  o  is  taken  to  be  zero  for  the  moment), 
although  the  mean  of  the  IG  distribution  is  m,  its  variance  is  the 
cube  of  p  divided  by  x . 

It  is  noted  that  the  a  parameter  is  a  redundant;  however,  it 

serves  as  a  location  parameter  which  allows  a  negative  range  of  vari¬ 
able  values.  It  may  be  useful  later  in  the  study  to  examine  the 

question  of  resolution  size.  The  behavior  of  p  and  x  by  will  be 

examined  looking  at  a  plot  of  x  versus  p  for  the  1x1  pixel  case. 

The  graphs  given  in  Figure  30  show  linear-linear,  log-linear,  and 
log-log  plots.  We  can  see  that  the  log-log  plots  separate  the  12 


scenes  somewhat  and  yet  significantly  groups  certain  fields.  At  the 
lower  left  portion  of  this  graph.  Clutter  1,  and  Clutter  2,  and  other 
grass  scenes  are  grouped.  These  scenes  are  all  quite  similar  in  that 
they  all  are  characterized  by  fairly  "smooth"  surfaces.  Next  come 
Fields  FI  and  F2  which  are  also  smooth,  yet  not  as  smooth  as  the 
grass  areas. 

Towards  the  upper  right  portion  of  the  graph  are  the  rougher 
scenes.  First,  the  Forests  1,  2,  and  3,  which  are  all  very  similar 
in  type,  mainly  deciduous  in  nature.  Fields  El,  E2,  61,  and  62  are 
all  of  very  rough  surfaces,  the  first  two  being  plowed  fields,  the 
latter  two  being  corn  fields. 

Thus,  using  the  16  characterization  parameters  u  and  x  suggests 
that  the  variation  from  lower  left  to  upper  right  represents  a  varia¬ 
tion  from  very  smooth  to  very  rough  scene  chartacteristics.  The 
linear-linear  and  the  log-log  plots  show  that  the  progression  from 
lower  left  to  upper  right  closely  approximates  a  straight  line  making 
a  45*  angle  with  the  axes.  At  the  moment  no  firm  explanation  for 
this  behavior  has  been  determined,  but  it  is  hypothesized  that  it 
may  be  connected  with  the  phenomena  of  surface  roughness  and/or  slope 
distribution.  So  it  is  very  possible  that  this  behavior  will  also 
depend  upon  the  resolution  question  mentioned  earlier  (that  is,  the 
grouping  of  pixel  sets). 

5.4  FUTURE  WORK 

The  results  presented  in  this  report  represent  a  sur-mary  of  the 
first  year  effort  of  the  analysis  of  SAR  generated  L-band  clutter 
data  (see  Figure  1  and  Table  1).  Efforts  during  the  next  year  of 
this  program  will  focus  on  the  analysis  of  this  data  as  described  in 
this  report.  In  addition,  emphasis  will  be  given  to  the  determina¬ 
tion  of  relationships  between  the  empirical  scattering  data  and  the 
measured  characteristics  of  the  scattering  surfaces.  An  expanded 
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data  set,  to  include  effects  of  incidence  angle  and  polarization, 
will  be  utilized  during  next  year's  effort. 

Additional  calibration  relationships  will  be  generated  using  SAR 
data  from  both  the  Long  Island  and  North  Carolina  reflector  arrays. 
Calibration  signal  generator  values,  which  are  recorded  on  each  data 
pass,  will  be  used  to  normalize  calibration  relationships  in  order 
to  improve  accuracy. 

Analysis  of  the  measured  data  describing  the  surface  properties 
will  also  be  continued  during  the  next  year.  The  ground  measurement 
data  that  has  not  been  considered  during  the  first  year  effort  will 
be  reduced.  Parameters  describing  the  slope,  radiis  of  curvative 
and  correlation  lengths  will  be  determined,  for  the  remaining  ground 
data.  Techniques  described  in  this  report  will  be  utilized,  in  addi¬ 
tion,  analysis  to  determine  improved  data  reduction  methods  will  be 
carried  out. 

Determinate  and  statistical  descriptions  of  the  ground  profile 
parameters  will  be  developed.  The  emphasis  will  be  directed  toward 
the  determination  of  surface  descriptors  which  correlate  with  the  16 
or  other  parameters  derived  to  discribe  the  scattering  data.  Consid¬ 
eration  will  be  given  to  the  size  of  the  data  set  and  resulting  sta¬ 
tistics.  Additional  data  samples  may  be  required.  The  deterministic 
and  statistical  parameters  required  for  use  in  the  several  candidate 
rough  surface  scattering  models  selected  will  be  calculated.  These 
models  will  be  used  to  determine  the  ability  to  describe  the  measured 
angular  scattering  response.  Also,  results  will  be  compared  with 
available  published  values  of  scattering  coefficients. 

SAR  data  from  the  test  sites  will  be  used  in  obtaining  additional 
parameter  values  for  the  16  distribution.  Values  obtained  for  the 
parameter  triplets  (i.e.,  u,  x,  a)  will  be  analyzed  to  determine 
distinct  characteristics,  properties  or  trends  as  a  function  of 
properties  of  the  scattering  surfaces.  Additional  data  sets  to  be 
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used  will  include  measured  values  over  a  range  of  incidence  angles 
for  a  number  of  the  test  sites.  The  first  analysis  will  be  confined 
to  horizontal  polarization  data.  After  this  initial  analysis  is 
completed,  vertical  polarization  data  from  the  same  test  sites  will 
be  analyzed  following  an  identical  procedure.  The  total  data  set 
will  then  be  analyzed  to  determine  characteristics  that  are  distinct 
between  horizontal  and  vertical  SAR  clutter  from  the  test  fields. 

There  are  two  separate  ideas  involved  in  the  estimation  of  the 
parameters  that  will  be  investigated  using  the  Inverse  Gaussian. 
First,  that  the  parameter  triplets  (a,  u,  x)  will  be  different  for 
distinct  terrain  types.  It  is  hypothesized  that  triplets  should  be 
distinct  for  distinct  incidence  and  polarization  views  of  the  same 
terrain  scene.  Second,  that  the  the  distribution  best  fitting  the 
data  is  the  IG  distribution.  The  first  can  be  accepted  or  rejected 
on  the  basis  of  plots  of  parameter  triplets  for  a  wide  variety  of 
terrain  data  which  is  available  at  ERIM. 

The  second  idea  is  more  difficult.  One  would  like  to  bring  the 
log-normal  distribution  into  play  in  the  examination  of  these  terrain 
data.  A  triplet  of  parameters  can  also  be  used  to  define  an  appro¬ 
priately  fitting  log-normal  distribution.  It  has  been  shown  for  ice 
data  that  in  many  cases  the  fitting  IG  and  log-normal  densities  can¬ 
not  be  visually  distinguished  from  each  other.  However,  in  that 
case,  the  IG  triplets  seemed  to  offer  wider  separation  for  different 
ice  types  than  the  log-normal  triplets.  To  make  a  final  choice  be¬ 
tween  the  IG  and  log-normal,  one  would  like  to  formulate  a  test  based 
on  a  ratio  of  maximum  liklihoods  for  these  two  distributions.  How¬ 
ever,  at  present  there  does  not  exist  a  table  of  critical  values  for 
such  a  test  involving  these  two  distributions.  And  to  establish  such 
a  table  will  require  some  sort  of  Monte  Carlo  procedure. 

It  is  planned  to  examine  goodness  of  fit  of  both  the  IG  and  log¬ 
normal  distributions  to  histograms  based  on  various  samples,  using 


the  Kolmogoroff-Smirnov  test.  The  analysis  and  tabular  data  neces¬ 
sary  to  carry  out  this  test  at  an  arbitrary  significance  level  al¬ 
ready  exist  in  the  statistical  literature. 

The  surface  height,  slope  and  radius  of  curvature  distributions 
will  be  used  in  a  statistical  analysis  similar  to  that  described 
above  for  the  scattering  data.  The  goal  of  this  part  of  the  analysis 
will  be  to  determine  correlation  between  surface  and  scattering 
statistics  within  the  limitation  of  the  data  sample  points  available. 
If  it  can  be  demonstrated  that  surface  profile  data  can  be  related 
to  scattering  characteristics  a  semi -empirical  model  may  result  for 
use  in  SAR  image  simulations  and  detector. 


6 

RECOMMENDATIONS 


The  present  program,  which  will  be  completed  at  the  end  of  FY 
1986,  was  intended  to  serve  as  an  initial  effort  in  collecting  and 
analyzing  calibrated  L-band  SAR  terrain  measurements.  This  effort 
was  limited  in  nature  because  only  a  small  set  of  SAR  data  were  util¬ 
ized.  This  limited  SAR  data  set  contained  a  sampling  of  the  various 
categories  of  terrain  types  which  may  be  considered.  Future  research 
efforts  in  this  area  should  be  focused  in  two  areas:  data  collection 
and  data  analysis. 

Future  data  collection  efforts  should  be  of  three  types:  (1) 
SAR  data;  (2)  scatterometer  data;  and  (3)  ground  measurements. 
Future  SAR  data  collections  would  utilize  the  ERIM  X-,  C-  and  L-band 
SAR  system.  Efforts  are  currently  underway  to  install  this  system 
into  a  U.S.  Navy  owned  and  operated  P-3  aircraft  (Shuchman,  et  al., 
1984).  As  was  accomplished  during  the  present  program,  SAR  data  of 
different  terrain  types  could  be  collected  on  a  target  of  opportun¬ 
ity  basis  during  test  flights  or  other  SAR  data  collection  programs. 

We  also  recommend  that  ground-based  L-band  scatterometer  measure¬ 
ments  of  various  terrain  types  be  collected.  Although  these  measure¬ 
ments  can  not  be  used  to  determine  the  statistical  distribution  of 
L-band  clutter  data,  this  type  of  data  would  be  useful  in  evaluating 
L-band  scattering  models.  This  data  type  is  much  more  easily  col¬ 
lected  than  airborne  SAR  data  and  can  be  used  to  help  verify  the 
scattering  models.  Scatterometer  data  would  be  collected  coincident 
with  the  SAR  data  collections  described  above  and  the  ground  measure¬ 
ments  discussed  below.  In  addition,  scatterometers  and  ground  mea¬ 
surements  could  be  collected  separately  from  the  SAR  missions. 

Finally,  a  ground  measurements  program  coincident  with  both  SAR 
and  scatterometer  data  collections  is  recommended.  The  ground  mea¬ 
surements  to  be  collected  include  dielectric  measurements  of  surface 
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materials,  soil  moisture,  soil  density  and  surface  roughness 
measurements. 

For  analysis  of  the  calibrated  L-band  terrain  measurements,  we 
recommend  that  research  be  conducted  in  three  separate  areas: 

1.  Further  evaluation  of  the  statistical  distribution  of  L-band 
terrain  data  and  the  correlation  to  ground  measurements; 

2.  Further  evaluation  of  L-band  clutter/scattering  models;  and 

3.  Evaluation  of  the  utility  of  L-band  SAR  image  simulation 
models. 

In  the  statistical  analysis  of  the  L-band  clutter  data,  efforts 
would  focus  on  the  extraction  of  radar  backscatter  signatures  from 
both  L-band  SAR  and  scatterometer  data  sets.  These  signatures  would 
then  be  correlated  with  the  surface  roughness  and  other  scattering 
statistics.  Finally,  we  could  continue  evaluation  of  the  statistical 
distributions  of  the  L-band  clutter  data. 

Using  the  L-band  terrain  measurements  and  surface  scattering 
statistics,  various  scattering  models  would  be  further  evaluated  for 
agreement  between  the  model  predictions  and  empirical  data. 

As  a  final  task,  we  recommend  that  the  utility  of  a  SAR  image 
simulation  model  being  developed  at  ERIM  be  evaluated.  An  airborne 
SAR  terrain  image  simulation  model  is  currently  being  developed  at 
ERIM  as  part  of  an  effort  being  sponsored  by  the  U.S.  Army.  Under 
the  recommended  future  effort,  we  would  concentrate  on  the  L-band 
portion  of  this  model.  The  model  would  be  upgraded  using  the  L-band 
terrain  measurements  obtained  during  the  recommended  program  as  well 
as  using  the  results  of  the  backscatter  model  evaluation.  The  imag¬ 
ing  model  would  be  evaluated  utilizing  SAR  imagery  available  at  ERIM. 
Finally,  a  point  target  simulation  model  which  has  been  developed  at 
ERIM  could  be  incorporated  into  the  terrain  model. 


APPENDIX 

EXAMPLES  OF  GROUND  MEASUREMENT  DATA 


This  appendix  presents  results  from  ground  measurements  made  dur¬ 
ing  the  North  Carolina  and  Long  Island  SAR  data  flights.  Figures  A-1 
to  A-12  includes  photographs  of  sites  in  the  Long  Island  area 
selected  for  initial  analysis.  Examples  of  measured  surface  height 
variation  after  digitization  of  the  raw  data  are  given  as  Figures 
A-13  through  A-24;  histograms  of  calibrated  a  values  from  example 
fields  are  given  in  Figures  A-25  through  A-30. 

Examples  of  results  obtained  from  reduction  of  the  surface  rough¬ 
ness  records  are  given  as  Figures  A-30  through  A-49.  These  include 
derived  examples  of  the  following  parameters 

Histograms  of  surface  slope  values  A-30  through  A-32 

Spectra  of  surface  slopes  A-33  through  A-37 

Spectra  of  surface  height  A-38  through  A-42 

Histograms  of  surface  height  A-43 

Histograms  of  surface  slope  A-44 

Histograms  of  radius  of  curvature  A-45  through  A-46 

Autocorrelation  of  surface  roughness  A-47  through  A-49 

Radar  data  from  grass  fields  in  the  reflector  test  array  sites, 
Peconic  River  Airport  and  from  fields  E,  F  and  G  have  been  used  in 
the  initial  statistical  analysis  reported  in  Section  5.  Examples  of 
ground  data  from  these  areas  are  included  in  this  appendix. 

In  addition  to  the  ground  data  included  in  this  appendix,  surface 
roughness  measuremets  are  available  from  other  fields  listed  in 
Table  7  for  the  North  Carolina  sites  and  Table  10  for  the  Long  Island 
sites.  Also,  photographs  showing  these  sites  are  available,  as  are 
tape  recordings  of  observations  made  during  the  7  September  1984  data 
collection  at  Peconic  River  Airport. 

Surface  roughness  measurements  were  made  at  the  Long  Island  test 
sites  during  the  period  27  August  to  6  September  1984.  All  soil 


samples  were  collected  on  7  September  1984  during  the  data  collection 
period  between  1200  to  1700  hours  EST.  Surface  roughness  measure¬ 
ments  were  made  of  the  North  Carolina  site  during  the  period  5  April 
to  10  April.  Also  soil  samples  were  obtained  on  8  and  9  April 
coincident  with  the  SAR  data  collection  times. 


(b) 

Figure  A-1.  a)  Ground  Photograph  of  Grass  Between  Reflector  Area  1 
and  2.  Peconic  River  Airport,  L.  I. 
b)  Photograph  of  48"  Square  Reflectors  -  Array  1 
Peconic  River  Airport,  L.  I. 


(b) 

Figure  A-2.  Ground  Photographs  of  Potato  Field  Near  Peconic  River 
Airport.  Designated  Field  3 

a)  Top  Photo  General  View  of  Field  Looking  West 

b)  Close  Up  Photograph  Showing  Row  Spacing  of  2' 10” 
and  Peak  to  Trough  Distance  of  8” 
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Figure  A-3.  Photographs  of  Two  Smooth  Fields  Northwest  of  Peconic 
River  Airport. 

a)  Smooth  Bare  Field  Designated  Field  N 

b)  Grass  Field  Designated  Field  A 


95 


✓ 


Figure  A-4.  Photographs  of  Field  C  Located  North  of  Peconic  River 
Airport  on  Sound  Avenue 

a)  View  Looking  South  From  Sound  Avenue 

b)  Close  Up  With  Lens  Cap  to  Show  Roughness  Scale 
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Photographs  of  Corn  Field  Designated  Field  G;  East  of  Field  F 

a)  View  of  Corn  Stalks  to  Show  Height,  Person  Shown  Is  6  Feet  Tall 

b)  View  of  Roughness  Board  in  Field  G 
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Figure  A-8. 


Photographs  of  Field  of  Cut  Corn  Designated  Field  H, 
East  of  Field  G 

a)  View  of  Field  H  Looking  Southeast 

b)  View  Showing  Roughness  Board 
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(b) 

Figure  A-9.  Photographs  of  Two  Grass  Fields  Located  North  of  Peconic 
River  Airport,  East  of  Landing  Road 

a)  Field  Designated  I,  Roughness  Board  Included  on  Photo 

b)  Field  Designated  K  Roughness  Board  Included  on  Photo 
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Figure  A-11.  Photographs  of  Smooth  Bare  Field  Located  Northeast  of 

Peconic  River  Airport  on  Middle  Country  Road.  Upper  Photo 
Shows  View  Looking  West  and  Lower  Photo  Shows  Close-up 
with  Roughness  Board.  Peak  to  Trough  2  inches.  Horizontal 
Peak  to  Peak  Distance  is  4  Inches.  This  Field  is 
Designated  P 


Figure  A* 12.  Photograph  Showing  Smooth  Grass  Field  Designated  Q, 
Located  West  of  Field  P.  Photo  Taken  Looking  North 
From  Middle  Country  Road 


SURFACE  HEIGHT  VERSUS  HORIZONTAL  DISTANCE 


Figure  A- 14.  Second  Roughness  Plot  of  Grass  Field  Site  2  Near  Reflector  Site 
Peconic  River  Airport 


Roughness  Plot  of  Grass  Near  Reflector 


FIKLD  ROUGHNKSS 

SURFACE  HEIGHT  VERSUS  HORIZON  l  AL  DIS  TANCE 


Figure  A- 16.  Field  ERT  From  Field  Designated  E,  Located  North  of  Peconic  River 
Airport  -  Roughness  Along  Peak  of  Furrows 


Designated  E,  Roughness  Normal  to  Furrows 


Figure  A- 18,  Roughness  Plot  From  Grass  Field  Designated  F,  North  of  Peconi 
River  Airport 


SURFACE  HEIGH 


Roughness  Plot  From  Field  H,  East  of  Field  E 


Smooth  Field  Designated 
Roughness  Measurement  T 


Roughness  Plot  From  Smooth 


Roughness  Plot  From  Smooth  Field  Designated  N,  East  of  Peconic 
River  Airport  on  Middle  Country  Road 
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Figure  A-25.  L-band  Histograms  of  Calibrated  Qq  Values  From  Test  Fields 


Figure  A-26.  L-band  Histograms  of  Calibrated  Oq '^*^*^®*  Fields 


Figure  A-27.  L-band  Histograms  of  Calibrated  Oq  Values  From  Test  Fields 


L-BAND  HISTOGRAM  -  FOREST  AREA  3 
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Figure  A-28.  L-band  Histograms  of  Calibrated  Oq  Values  From  Test  Fields 


A-29.  L-band  Histograms  of  Calibrated  Values  From  Test  Fields 


Figure  A-30.  Histogram  of  Surface  Slope  Field  EG  -  Potatoe  Field  E,  North 
of  Peconic  River  Airport 


Figure  A-33.  Field  H  -  Cut  Corn,  North  of  Peconic  River  Airport 


North  of  Peconic  River  Airport 


Reflector  Field  RFS3R  Grass  Field  Near  Reflector  Site  3,  Peconic 
River  Airport 


Reflector  Field  RFS2  Grass  Field  Near  Reflector  Site  2,  Peconic 
River  Airport 


ield  North  of  Peconic  River  Airport 


Figure  A-36.  Field  E  -  Potato  Field  North  of  Peconic  River  Airport 
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Field  E  -  Potato  Field  North  of  Peconic  River  Airport,  Along 
Peak  of  Furrows 


Figure  A-40.  Field  F  -  Grass  Field  North  of  Peconic  River  Airport 


Figure  A-41.  Field  H  -  Cut  Corn,  North  of  Peconic  River  Airport  -  E-W  Line 


Figure  A-42.  Field  H  -  Cut  Com,  North  of  Peconic  River  Airport  -  N-S  Line 
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Figure  A-43.  Histogram  of  Surface  Height  -  Grass  Field  F,  Height  Values  Given  In  Inches 


Figure  A-44.  Histogram  of  Surface  Slope  -  Field  E,  Along  Peak  of  Furrows 


Figure  A-45.  Histogram  of  Surface  Radius  of  Curvature,  Field  E,  Along  Peak  of 
Furrows  -  Radius  Values  Given  In  Inches 


Figure  A-49.  Potato  Field  E  -  Furrows 
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